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INTRODUCTION 


The  finite  element  code  may  be  used  to  analyze  plane  strain  and 
axisymmetric,  quasi-static,  soil  problems,  including  consolidation  effects.  The 
soil  may  be  modeled  using  either  linear  elasticity  or  the  "bounding  surface 
plasticity  model  for  cohesive  soil".  The  program  is  written  in  modular  form  so 
that  other  soil  models  can  be  easily  incorporated.  The  theory  underlying  the 
analysis  is  described  fully  in  the  accompanying  report  [1]. 

The  user's  manual  is  divided  into  six  parts:  Introduction,  Input,  Output, 
Explanatory  Comments,  References  and  Example.  The  Input  section  gives  in 
outline  form,  the  sequence  of  information  required  to  describe  the  problem  to 
be  analyzed.  Only  the  briefest  notes  of  explanation  are  included  in  this  section. 
Until  an  analyst  becomes  familiar  with  the  program  he  or  she  will  need  to  refer 
to  the  last  section  of  the  manual  where  detailed  explanations  and  examples  are 
given.  The  manual  assumes  general  familiarity  with  finite  element  methods; 
novices  in  this  area  are  referred  to  a  standard  text  such  as  [2]. 
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n.  INPUT 


The  required  input  data  is  entered  by  means  of  the  following  sequence 
of  records  (cards): 

At.  Title  Card  (!SA»)t 

Any  information  that  is  to  be  printed  as  the  title  of  the  problem. 
A2.  Control  Card  <L5,  F5.0,  15,  F5.0,  15,  F5.0,  715) 


Columns 

5 

iquit 

(T(rue)  -  stop  i  analysis  after 

[p(alse)  -  do  not  stop)  roesh  S'nera,i,ln 

6-10 

GRIDV 

= 

grid  generation  parameter  (0.0  -"isopara¬ 
metric"  grid,  1.0  -  "Laplacian"  grid) 

15 

IFLOW 

10  unsaturated  problem  -  no  water  flow 
|1  saturated  problem  -  water  flow 

16-20 

ei 

= 

parameter  controlling  numerical  inte¬ 
gration  in  time 

25 

MTYPE 

s  < 

[0  -  plane  stress  analysis 
,  1  -  plane  strain  analysis 
[2  -  axisymmetric  analysis 

26-30 

a 

s 

parameter  controlling  "reduced"  integra¬ 
tion  of  volume  term 

Upper  bounds  on  dimensions  (used  to  establish  dynamic  storage 
allocation,  all  values  except  NCOFMX  >  0): 


31  -  35 

NFUNMX 

> 

no.  of  history  function  specifications. 
Section  B1 

* 

1 

s 

IFUNSZ 

> 

greatest  no.  (M)  of  points  used  to  describe 
a  given  history  function,  Section  B1 

41  -  45 

MATMX 

> 

no.  of  materials.  Section  B2 

46-50 

NCOFMX 

> 

no.  of  initial  state  descriptions.  Section 
B3 

51  -  55 

NPTMX 

> 

the  largest  node  number.  Section  B4 

56-60 

NELMX 

> 

no.  of  elements.  Section  B5 

61  -  65 

NDSPMX 

> 

no.  of  node  point  specifications.  Sec- 
tion  B6 

Columns 


A3.  Gi 


1  -  10 

S 

=  acceleration  of  gravity  (magnitude) 

11  -  15 

% 

=  history  function  number  for  magnitude  of 
gravity 

16  -  25 

=  angle  (in  degrees)  that  gravity  makes  with 
negative  y(z)  axis 

26-30 

=  history  function  number  for  9 

g 

A#.  Nonlinear  Analysis  Card  (13.  F5.0.  313,  2E10.3): 

Specification  of  desired  iteration  options  (leave  blank  for  a  nonlinear 
problem  ) 


Columns 

3  NONUN 


!0  linear  problem 
^[nonlinear  problem-job 


terminated 

is  [if  convergence 
is  not!  does  not 
,  occur 


6-10  0.GH1<0.5  =  parameter  controlling  Newton-Raphson  approx¬ 

imation  (0.0  gives  the  tangent  stiffness 
method;  0.5  gives  the  method  of  successive 
approximation) 


11  -  15  UMAX 

16-20  IRPET 

25  ITFAC 

26  -  35  FL 

36-45  ERMAX 


maximum  number  of  iterations  permitted  in 
any  single  solution  increment  (default 
value  =  5) 


5f 


reform  stiffness  matrix  every 


i  iteration 
K-th  iteration 


fT(rue)  \ 
F(alse)| 


variable  acceleration  factor  applied  to 
solution  vector  components 


* 


places  limits  of  1/FL  (  )  >.  FL  on  the 
acceleration  factor  when  ITFAC  =  T  (default 
value  =  0.3) 

convergence  criterion  for  the  solution  vector 
(default,  value  =  0.01) 


after  the  2nd 


3 


A  card  with  a  1  punched  in  column  1  followed  by  Of  no  history  function 
specifications  are  required,  this  section  is  omitted  entirely): 

History  Function  Descriptions:  The  following  cards  are  required  for  each 

distinct  function  (History  functions 
numbered  -3,  -2,  -I  and  0  are  explicitly 
defined  in  the  program,  see  Explanatory 
Notes,  and  thus  no  input  is  required): 

1st  Card  (IX.  I».  j Sk 


Columns 


2-5  IH 
6-10  M 
2nd  CardU)  (8E10.3): 


function  number  (>0) 

number  of  points  needed  to  define  the  function 


As  many  cards  as  needed  to  specify  the  M  pairs  of  values  (F  ,  t  ) 
which  defines  the  function.  The  initial  card  contains  the  values 
*1»  F2».  t2,'lM  F4»  t4*  Subsequent  cards,  if  required  (M  >  4),  contain 
tne  values  F5,  tj,...,  F,^,  t^. 


jt-WTAt'  “.Wl? 


A  card  with  a  2  punched  in  column  1,  followed  by: 


p  The  following  information  must  be  supplied 
for  each  distinct  material: 


'  & 

:■$ 


& 

18 

1  d 


<$• 

M 


|-:S 

isa 


-V* 

?  * 


'i\ 


\i  1 

\  i 


"t 

*  V 


lit  Card  (IX.  14.  15.  6E10.3): 
Columns 

2-5  NMAT 

10  ITYP 


=  material  number 


!1  -  isotropic  linear-elastic 
2  -  anisotropic  linear-elastic 
3  -  bounding  surface  plasticity  model  for 
cohesive  soil 


11 

-  20 

ps 

s  soil  density+ 

21 

-  30 

pf 

s  fluid  density ++ 

31 

-  40 

r 

=  bulk  modulus  for  fluid  and  soil  particles 

41 

-  50 

le* 

*11  . 

1 

51 

-  60 

•ft 

[  =  effective  soil  permeability  coefficients 

61 

-  70 

•5* 

2nd  Card  (8E10.3); 


Columns 

ITYP  =  1 

ITYP  =  2 

ITYP  =  3 

1  -  10 

G 

Dn 

X 

11-20 

V 

DJ2 

K 

21-30 

°13 

Mc 

31  -  40 

DI* 

Rc 

41-50 

°22 

Ac 

51-60 

D23 

T 

61  -  70 

°24 

P1 

71-80 

°33 

v  or  G 

♦  If  the  acceleration  of  gravity  (g  ec..  xi  A3)  is  taken  as  unity  then  p 
and  Pg  are  unit  weights. 

♦+  If  it  is  desired  to  use  Excess"  not  total  pore  water  pressure  then  p,  is 
set  equal  to  zero  -  see  explanatory  comments. 


C’4 


3rd  Card  (8E.10.3)  -  required  only  if  ITYP  >  1 


Columns 

ITYP  =  1 

ITYP  =  2 

ITYP  =  3 

1  -  10 

°34 

Pa 

11-20 

°44 

r  (Duplicates  value  on 

21  - 

30 

m 

31  - 

40 

h 

41  - 

50 

c 

h2 

51  - 

60 

n  =  Me^Mc 

61  - 

70 

»  -  V>c 

71  - 

80 

r  *  V’c 

4th  Card  (3E10.3)  -  required  only  if  ITYP  =  3: 

Columns  ITYP  =  1  ITYP  =  2 

ITYP  =  3 

1  - 

10 

a  =  Ae/Ac 

11  - 

20 

C 

"1st  can 


A 


B3.  A  card  with  a  3  punched  in  column  1,  followed  by: 


Initial  State  Information:  The  following  information  must  be  supplied  for 

each  non-trivial  initial  state  (this  section  is 
omitted  if  no  information  is  required)* 


1st  Card  (IX.  14,  6E10.3); 
Columns 


2-3 

ISNO 

= 

state  number 

6-15 

al) 

- 

initial  effective,  vertical  stress  distr- 

16  -  25 

41 

ibution,  in  the  form  ay  =  a^  +  ajy 

26-35 

b>l 

s 

initial  effective,  horizontal  stress  distri- 

36  -  45 

bution,  in  the  form  =  bj  +  b2y 

46  -  55 

c,i 

s 

initial  pore  water  pressure  distribution, 
in  the  form  h  =  Cj  ♦  c^f 

56  -  65 

c2 ; 

aid  Card  (4E10.3) 

1-10  dj  ) 

11-20  d2) 

21  -  30  ejJ 

31  -  40  e2) 


initial  void  ratio  distribution,  in  the  form 

ej  s  dj  +  djy 

initial  preconsolidation  pressure  distri¬ 
bution,  in  the  form  PQ  =  e^  + 


*  The  initial  state  of  o  =  o.  =  h  =  e.  =  P  =  0,  is  given  the  number  zero 
and  is  built  into  the  program.  The  equations  a,  +  a2y >  etc.  assume  plane 
conditions.  For  axisymmetry,  the  vertical  coordinate  direction  is  z  instead 
of  y,  and  the  distributions  are  then  of  the  form  oy  =  a^  ♦  a^z,  etc. 


A  card  with  a  %  punched  in  column  1,  followed  by: 

Node  Point  Array  (IX,  14,  2E10.3,  15,  3E10.3):  As  many  cards  as  are 

necessary  to  specify  the 
locations  of  all  nodes  in 
the  system  which  are  not 
to  be  generated  by 
means  of  the  "interior" 
generation  scheme: 

Columns 


NJ 

i 


2  -  5 

N 

= 

node  point  number 

6-15 

X 

s 

x(r)  -  coordinate 

16  -  25 

Y 

s 

y(z)  -  coordinate 

26-30 

INC 

s 

numbering  increment 

o 

1 

D 

= 

spacing  ratio 

41  -  50 

XC  i 

e 

coordinates  of  some 
point  on  the 

51-60 

YC  ) 

interior  of  the 

circular  arc 


quantities 
associated 
with  the 
straight  and 
curved  line 
generation 
options 


1 


i;j 


1 


.  «  -  .’n  .  ’•  *  \ 


B5.  A  card  with  a  5  pinched  In  column  l,  followed  by: 


Element  Array  (IX.  14.  flSk  As  many  cards  as  are  necessary  to  specify 

all  elements  in  the  system:* 


Columns 

2-5 
6-10 
11  -  15 
16  -  20 


•ry 

IN 

1 

•■w 

fM 

MN 

26-30 

BNO 

31  -  35 

NMISP 

• 

o 

INCRP 

41  -  45 

NMB 

46-50 

INCR 

the  numbers  of  the  four  node  points 
which  describe  the  quadrilateral  or 
triangular**  element  (reading  counter¬ 
clockwise  around  the  element) 

material  number  (corresponding  to  the 
appropriate  material  description  of 
section  B2) 


initial  state  number  (corresponding  to  the 
appropriate  initial  state  description  of 
section  B3) 


number  of  additional 
elements  in  the  layer 

numbering  increment 
for  elements  within 
the  layer 


number  of 
additional  layers 


quantities 

associated 

with  the 

element 

generation 

option 


numbering  increment 
for  the  layers 


*  The  order  of  the  element  cards  need  bear  no  relation  to  the  actual  location 
of  the  elements  within  the  body.  The  order  will  determine  the  assigned 
"element  numbers." 


For  a  triangular  element  the  forth  node  number  is  set  equal  to  the  first. 


B6.  A  card  with  a  6  punched  In  column  1,  louowed  oys 


Node  Point  Specification  Array  (IX.  14.  303.  12.  E1Q.3).  El 0.3.  213.  2F5Ak 


As  many  cards  as  are  necessary  to  specify  known 
nodal  displacements,  loads,  water  flows  and  pore  water 
pressures. 


Columns 


2-5 

N 

6-8 

!Hi 

10 

*1 

11  -  20 

V1 

21  -  23 

"*2 

25 

*2 

26  -  35 

V2 

36-38 

m, 

40 

*3 

41-50 

V3 

51  -  60 

®n 

61  -  65 

N 

66-70 

INC 

71  -  75 

PN 

76  -  80 

PN* 

Node  point  number 


history  function  number  (Section  Bl)  for 
the  l-coordinate  direction 

“{indicates  that  a  known  j,<S"Ketnent  j 
is  specified  in  the  1-coordinate  direction 

magnitude*  of  the  specified  jdf^p"cem„t 
for  the  1 -coordinate  direction 

history  function  number  (Section  Bl)  for 
the  2-coordinate  direction 

“{indicates  that  a  known  |f"^acelT,ent| 
is  specified  in  the  2-coordinate  direction 

magnitude*  of  the  specified  j^ceitwit 
for  the  2-coordinate  direction 


history  function  number 

^vindicates  that  a  known 
is  specified 


water  flow 
pore  water  pressure 


magnitude*  of  the  specified 


water  flow 
pore  water  pressure 


angle  (in  degrees)  between  the 
Xj-axis  and  the  x(r)-axis 


final  node  point 
in  the  sequence 

numbering  increment  for  node 
points  in  sequence 

values  of  linear  varying 
pressures  applied  at  points 
N  and  N'  respectively** 


quantities 

associated 

with  the 

node  point 

specification 

generation 

option 


*  In  all  cases  the  actual  value  of  the  prescribed  quantity  is  the  product 
of  the  "magnitude"  and  the  value  of  the  specified  "history  function" 

**  The  history  function  number  of  P  must  be  prescribed  in  IH,  and  IHj 


B7.  A  card  with  a  7  punched  in  column  1,  followed  by: 


Solution  History  Segment  Information  (IX.  II,  2E10.3):  One  card  for 

each  history 
segment  into 
which  the  incre¬ 
mental  analysis 
is  subdivided:* 


Columns 

2-5 

NMB 

number  of  solution  (time)  increments  into 
which  the  history  segment  is  subdivided 

6-15 

TIME  = 

time  at  the  end  of  the  history  segment 

16  -  25 

D 

incrementing  ratio  controlling  the  time- 
step  lengths  within  the  history  segment 
(default  value  =  1.0) 

End  Card  Ql)t 

A  card  with  an  8  punched  in  column  1  to  denote  the  end 
of  the  input  data  for  the  given  problem. 

The  above  sequence  of  cards  A1  ♦  C  are  repeated  for  each  additional 
analysis  in  the  "stack". 


*  note  that  the  analysis  begins  at  time  t  =  0 


m.  OUTPUT 


The  output  from  the  program  consists  of  am  echo  print  of  material 
properties  and  solution  parameters,  the  generated  node  and  element  data, 
messages  for  detected  data  errors,  and  finally  for  each  time  step  the  problem 
solution.  When  data  errors  are  detected,  the  program  aborts  the  job  after  the 
printing  of  the  input  data  and  proceeds  to  the  next  job  in  the  stack  of  data. 
The  message  "data  error  in  element  I"  is  triggered  by  a  computed  negative  area 
for  one  or  more  of  the  four  triangles  formed  by  three  successive  vertices  of 
the  element.  Such  a  negative  area  is  a  result  of  either  an  improperly  shaped 
element  resulting  from  a  poor  mesh  selection,  a  data  error  in  entering  node 
point  coordinates  (often  for  one  of  the  nodes  defining  the  element  in  question) 
or  element  data  (including  entering  the  nodes  CW  instead  of  CCW  -  assuming 
the  coordinate  system  to  be  right-handed). 

The  print  out  of  the  node  point  specifications  includes  any  concentrated 
node  point  forces  (in  x-y  components)  resulting  from  specified  surface  pressures. 

The  printed  values  of  strains,  stresses,  etc.,  at  a  given  time  step,  are 
the  values  accumulated  to  that  point  in  time  including  initial  values.  The 
stresses  are  effective  stresses  (tension  positive).  The  pore  water  pressure  (units 
of  stress  -  compression  positive)  will  be  either  total  pressure  or  "excess"  pressure 
depending  on  user  preference,  see  Section  B2  in  part  IV. 

The  headings  for  the  solution  output,  are  self  explanatory  with  the  possible 
exception  of  h  which  denotes  the  pore  water  pressure. 


IV.  EXPLANATORY  NOTES  REGARDING  THE  INPUT 


General  Comments: 

The  following  comments  hold  for  Sections  B1  B7.  For  example,  consider 
Section  B4.  The  first  card  in  the  Section  contains  the  number  4  punched  in 
column  one;  it  has  no  other  information  (only  one  such  card  is  used).  It  is 
followed  by  as  many  additional  cards  as  are  needed  to  define  the  node  point 
data.  The  first  column  of  each  of  these  cards  is  left  blank. 

It  is  the  responsibility  of  the  user  to  maintain  consistent  units.  The  units 
used  to  describe  gravity  (Section  A3),  and  the  material  properties  (Section  B2) 
must  be  consistent  with  those  used  to  describe  the  initial  state  (Section  B3), 
the  geometry  of  the  body  (Section  B4),  and  the  node  point  specifications  (Section 
B6).  The  solution  is  expressed  in  the  same  units  as  the  input. 

Because  the  bandwidth  NBAND  of  the  simultaneous  equations  is  determined 
by  the  numbering  of  the  nodes,  an  optimal  node  numbering  scheme  is  required 
to  minimize  the  computational  cost  of  a  given  finite  element  analysis.  The 
bandwidth  resulting  from  a  given  numbering  scheme  is  computed  in  the  following 
manner: 

i)  Denote  the  span  for  any  two  nodes  of  a  given  element  as  N.,  where 
N.  is  equal  to  the  absolute  difference  in  the  node  numbers. 

ii)  Denote  the  maximum  value  of  Nj  for  a  given  element  j  as  NE.. 

iii)  Considering  all  elements  in  the  system,  denote  the  maximum  value 

of  NE.  as  NE  . 
j  max 

iv)  The  bandwidth  is  then  given  by  the  expression  NBAND  a  (2  + 
IFLOW)*(NEmax  +1) 

Since  NEmax  js  directly  related  to  the  bandwidth  of  the  simultaneous 
equations,  in  numbering  the  nodes  it  is  this  quantity  that  should  be  minimized. 
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Section-by-Section  Comments: 

The  section  numbers  used  below  correspond  to  the  section  numbers  of 
part  D.  INPUT,  thus,  in  order  to  find  information  concerning  the  input  for  B6 
(Node  Point  Specification  Array)  the  reader  should  refer  to  Section  B6  below. 
In  addition,  within  a  given  section  items  called  out  in  the  input  are  typed  in 
bold.  For  example  input  items  6  and  ITMAX  which  are  required  for  input  Section 
A4  are  type  in  bold  where  they  are  discussed  below  in  Section  A4.  The  theory 
underlying  the  analysis  is  only  superficially  treated  here;  for  a  more  complete 
discussion  the  reader  is  referred  to  [1]. 

Al.  Title  Card 

The  title  serves  to  identify  the  particular  problem  under  consideration. 

A2.  Control  Card 

If  a  T(rue)  value  is  specified  for  IQUIT  the  analysis  terminates  after  the 
mesh  has  been  generated  and  printed.  This  option  should  be  used  for  the  first 
run  of  a  large  problem  in  order  to  avoid  wasting  computer  time  analyzing 
incorrect  data.  If  data  for  several  problems  is  contained  in  the  stack,  the 
program  skips  the  time  history  data  for  the  terminated  job  and  proceeds  to  the 
next  problem. 

For  the  precise  meaning  of  the  grid  generation  parameter  GRIDW  the 
reader  is  referred  to  13]  (GRIDW  =  1.0  -  w,  where  w  is  defined  in  13]).  In 
general  a  value  of  GRIDW  =  0.0  is  recommended;  for  those  very  rare  cases 
where  this  results  in  a  singular  set  of  equations  for  the  grid  generation  process, 
a  value  of  .05  is  recommended. 

The  code  IFLOW  distinguishes  between  saturated  conditions  where  water 
flow  occurs  (or  a  potential  for  water  flow  exists  -  ideal  undrained  conditions) 
and  unsaturated  conditions. 
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When  IFLOW  =  1  the  soil  density  (or  unit  weight  if  the  acceleration  of 
gravity  is  taken  to  be  unity)  pg  specified  in  Section  B2,  refers  only  to  the  soil 
skeleton  (unsaturated  soil).  The  printed  stresses  are  the  "effective  stresses?'  and 
must  be  supplemented  by  the  pore  water  pressure  to  obtain  the  total  stresses. 
If  it  is  desired  to  exactly  model  "ideal  undrained  conditions?1  (no  movement  of 
water),  the  effective  permeability  of  the  soil  should  be  set  equal  to  zero  (Section 
B2). 

When  IFLOW  =  0  the  soil  density  must  include  the  mass  (or  weight) 
of  any  water  present  in  a  partially  saturated  soil  (the  pore  water  pressure  is 
assumed  to  be  zero  and  water  is  assumed  not  to  flow).  The  printed  stresses 
are  total  stresses. 

Conditions  where  part  of  the  soil  mass  is  unsaturated  and  part  is  saturated 
can  be  modeled  by  specifying  for  the  unsaturated  soil  a  very  small  bulk  modulus 
T  for  the  water  (and  soil  particles  -  Section  B2). 

The  parameter  0  j  determines  the  approximation  used  for  the  time 
derivatives  in  the  governing  equations  (see  [1]);  values  between  0.5  (Crank- 
Nioolson)  and  .67  (Galerkin)  are  recommended  [2], 

For  axisymmetric  analyses  (MTYPE  =  2)  all  x  and  y  notations  refer  to  r 

and  z. 

The  parameter  a  determines  the  finite  element  approximation  used  for 
measuring  volume  change  (see  [1]).  When  IFLOW  =  1,  a  value  of  0.0  is 
recommended  unless  solution  oscillation  is  a  problem  in  which  case  a  value  of 
.1  is  useful.  Except  for  nearly  incompressible  linear  elastic  materials,  when 
IFLOW  =  0  a  value  of  1.0  is  usually  preferable. 

All  arrays  in  the  program  whose  dimensions  are  problem  dependent,  are 
dynamically  dimensioned.  The  values  MATMX, .  .  .  NDSPMX  contain  information 
for  this  purpose,  all  these  quantities,  with  the  exception  of  NCOFMX,  must  be 


greater  than  zero.  The  value  of  MATMX,  etc.  are  upper  bounds  and  thus,  unless 
it  is  desired  to  absolutely  minimize  storage  requirements,  need  not  be  equal  to 
the  actual  number  of  specified  materials,  etc.  When  specifying  the  values 
NELMX  and  NDSPMX,  it  must  be  remembered  to  count  those  elements  and 
specifications  which  are  included  by  means  of  the  generation  options. 

In  the  dynamic  dimensioning  of  the  program,  separate  arrays  were  used 
far  interger  and  floating  point  numbers  in  order  to  avoid  difficulties  for  computers 
that  use  different  word  lengths  for  the  two.  The  program  has  been  coded  so 
that  16  bit  integers  may  be  used  if  desired. 

The  dimensions  of  the  program  are  controlled  by  two  quantities  "long" 
and  "longi"  specified  in  "PARAMETER"  statements  at  the  beginning  of  the 
program;  these  quantities  must  satisfy  the  following  inequalities, 
longi  _>  NFUNMX  +  NPTMX  +  5*NELMX  +  (3+IFLOW)*NDSPMX  +  1 
long  _>  26*MATMX  +  2*NFUNMX*IFUNSZ  ♦  NCOEF*(NCOFMX  +  1) 

♦  [2  ♦  4*(2  ♦  DFLOW)l*NPTMX  +  (3  +  IFLOW)*NDSPMX 
♦  LONGEQ 

Where  LONGEQ  is  the  space  set  aside  for  solving  the  system  of  equations  by 
means  of  a  block,  constant  bandwidth  equation  solver.  If  the  bandwidth  of  the 
equations  is  denoted  as  NBAND,  then  the  minimum  value  for  LONGEQ  is 
NBAND*NBAND  (only  a  single  equation  would  be  contained  in  each  equation 
block);  if  it  is  desired  to  solve  the  equations  entirely  in  core  then  it  must  have 
a  value  NBAND*(2  ♦  IFLOW)*NPT.  In  general  it  is  recommended  that  LONGEQ 
exceed  the  minimum  by  at  least  30%.  The  calculation  of  the  bandwidth  NBAND 
is  discussed  in  the  general  comments  at  the  beginning  of  this  part  of  the  manual. 
A3.  Gravity  Card 

Gravity  g  can  be  input  either  in  terms  of  the  acceleration  units  appropriate 
to  the  system  of  units  selected  for  the  problem  (32.2  ft/sec  for  English  units) 
or  in  terms  of  multiples  of  the  acceleration  of  gravity  at  sea  level  (i.e.  g  =  1 
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for  a  field  structure);  the  corresponding  meanings  of  and  (Section  B3) 


would  be  mass  densities  in  the  first  case  and  unit  weights  in  the  second.  That 


is  the  product  pg  must  have  units  of  weight  per  unit  of  volume. 


The  histories  of  the  magnitude  g  and  direction  0  of  gravity  are  specified 


by  the  history  function  numbers  (Section  Bl)  IH  and  IH*.  A  pre-existing  gravity 

O 

loading  of  a  field  deposit  can  be  modeled  by  initializing  the  stresses  and  pore 


water  pressure  (Section  B3)  to  their  proper  values  and  setting  IH  =  IHq  =  -3. 


The  history  of  the  effective  gravity  loading  on  a  centrifuge  model  during  "spin-up" 


can  be  modeled  by  describing  in  Section  Bl  a  history  function  corresponding  to 


the  centrifuge  velocity  history  for  the  test;  in  the  case  of  a  fixed  bucket  both 


g  and  0  would  vary  with  time,  while  for  a  swing -up  bucket  only  g  would  vary. 

O 


M.  Nonlinear  Analysis  Card 


For  a  linear  elastic  problem,  NONLIN  is  set  equal  to  zero  and  the  rest 


of  the  card  is  left  blank.  For  problems  using  the  bounding  surface  plasticity 


model,  NONLIN  is  set  equal  to  1  or  2  depending  on  whether  the  analysis  should 


be  terminated  or  not,  if  convergence  is  not  achieved  in  a  given  time  step. 


The  factor  0  determines  the  approximation  to  be  used  for  the  Jacobian 


in  the  Newton-Raph son's  iteration  for  the  nonlinear  problem,  for  details  the 


reader  is  referred  to  [4,  5],  It  is  expected  that  a  value  of  0.0  will  in  most 


cases  give  the  best  results. 


The  frequency  of  updating  the  stiffness  matrix  during  the  iteration  process 


is  controlled  by  the  value  of  IRPET  [4];  for  initial  uses  of  the  program  a  value 


of  zero  would  appear  to  be  appropriate.  The  values  of  ITFAC  and  FL  control 


the  use  of  acceleration  factors  applied  to  the  components  of  the  solution  vector 


[41;  for  initial  use  of  the  program  it  is  suggested  that  ITFAC  =  F.  Finally, 


the  default  value  of  .01  for  the  convergence  limit  ERMAX  would  appear  to  be 


adequate  for  most  problems. 


Bl.  History  Ruction  Descriptions 

The  time  dependence  of  all  input  quantities  (i.e.,  magnitude  and  direction 
of  gravity,  and  specified  node  point  displacements,  flows,  loads  and  pore  water 
pressures)  are  specified  by  means  of  appropriate  'history  functions".  The  program 
has  built-in  four  such  functions  numbered  -3  -*•  0,  i.e., 


i)  IH  =  -3 


U)  IH  =  -2 


Specifies  a  in  it  value  and  a  zero  incremental  value 
for  aU^  times,  Figure  la. 

Specifies  a  zero  value  and  a  zero  incremental  value 


iv)  IH  =  0: 


for  all  times  -  Figure  lb. 

iii)  IH  =  -1  Specifies  all  incremental  values  equal  to  1.0.  The 

incremental  values  are  taken  to  be  equal  regardless 
of  the  relative  lengths  of  the  time  steps  specified  in 
Section  B7.  The  resulting  history  functions  for  the 
cases  of  equal  and  variable  length  time  steps  are 
illustrated  in  Figure  lc. 

iv)  IH  =  0:  Specifies  a  step-function  at  time  t  =  0;  that  is  a 

quantity  using  this  function  is  applied  entirely  Airing 
the  first  solution  increment,  Figure  Id. 

In  addition,  the  user  may  described,  by  means  of  the  input  to  Section  Bl, 
as  many  more  history  functions  (numbered  1-0  as  needed;  an  example  of  such 
a  function  is  given  in  Figure  2. 

For  a  particular  history  function,  linear  interpolation  is  used  to  identify 
die  AF  which  corresponds  to  a  given  time  increment  At.  For  solution  times 
beyond  the  last  specified  point  V  the  final  history  segment  is  extended 
indefinitely. 

When  a  magnitude  V  and  a  history  function  number  IH  are  specified  in 
Section  B6  (or  A3)  for  some  given  external  agent,  then  in  the  solution  interval 


IS 


At  an  incremental  value  of  the  quantity  equal  to  V*AF  is  applied,  where  AF 
corresponds  to  history  function  IH. 

B2.  Material  Properties  Array 


The  units  of  the  material  properties  must  be  consistent  with  the  units 
used  to  describe  the  geometry  of  the  body  and  the  magnitudes  of  the  applied 
loads. 

The  material  number  NMAT  serves  as  an  identifier  for  use  in  Section  B5 
to  assign  a  particular  material  description  to  a  group  of  elements.  In  the  current 
version  of  the  program  three  types  of  material  descriptions  are  permitted,  i.e., 
isotropic  or  anisotropic,  linear  elastic  and  the  bounding  surface  plasticity  model 
for  cohesive  soils.  Additional  material  models  cam  be  easily  added  to  subroutine 
PROPTY  by  extending  the  two  key  "Block  IP*  statements  as  indicated  in  the 
program  by  comment  statements. 

As  noted  previously  (in  Section  A3),  the  units  of  ps  and  must  be 
compatible  with  the  units  selected  for  gravity  "g". 

Flow  problems  (IFL0W=1)  can  be  stated  either  in  terms  of  total  or  excess 
pore  water  pressure.  In  the  first  case  Pj  must  be  set  equal  to  the  fluid  density 
(or  unit  weight  -  see  previous  paragraph);  in  the  second  case  it  is  set  equal  to 
zero. 

The  quantity  T  can  be  viewed  either  as  the  combined  bulk  modulus  of 
the  soil  particles  and  the  pore  water,  or  as  a  penalty  number  imposing  an 
assumed  incompressibility  condition  for  these  components  [1,  2,  6J.  In  the 
absence  of  experimental  evidence,  the  bulk  modulus  for  water  (3.2  x  105  psi, 
2.2  x  10^  N/m^)  may  be  used  for  T. 

The  "effective"  permeability  coefficients  kjj  appear  in  Darcy’s  law^  when 

Uftits  of  pressure  (not  head)  are  used  for  the  pore  water  pressure;  their 

'  In  terms  of  excess  pore  water  pressure  it  has  the  form 
„*  3  h  , *  a  h  *  _ 


relationships  to  the  permeability  coefficients  commonly  used  by  civil  engineers 
and  those  used  by  physicists  are  discussed  in  [1]. 

For  isotropic,  linear  elasticity  E  and  \>  denote  Young's  modulus  and  Poisson's 
ratio  respectively.  The  linear,  anisotropic  elastic  law  is  written  in  the  form: 


T 

For  plane  (x-y)  conditions  the  stress  vectors  is  <o'  o'  o'  o'  >  ,  while 

x  y  z  xy 

T 

for  axisymmetry  it  is  <o',  o *,  a g,  t  '^>  . 

The  meanings  of  the  several  parameters  describing  the  bounding  surface 
model  are  described  in  detail  in  ref  17-11];  values  for  particular  soils  may  be 
found  in  [8,  11];  a  summary  of  information  is  given  in  Table  1.  In  order  to 
keep  the  input  for  the  model  exactly  as  described  in  [10],  the  parameter  T  is 
retained  even  though  it  is  a  duplication  of  previous  input  (the  second  value  is 
not  used). 

B3.  Initial  State  Information 

The  information  in  this  section  is  used  to  establish  the  initial  state  of 

the  soiL  The  values  of  h  specified  in  this  section  are  used  directly  to  initialize 

the  pore  water  pressure  in  the  elements  and  indirectly  to  initialize  it  for  the 

nodes;  see  Section  B5.  It  is  extremely  important  to  note  that  ah  and  oy  are 

"effective"  stresses  (total  stress  minus  pore  water  pressure).  It  is  assumed  that 

the  initial  horizontal  stress  is  the  same  in  all  directions  and  the  initial  shear 

stress  (t  or  t  )  is  zero.  For  linear  elasticity  problems  the  initial  stresses 
xy  rz 

may  be  taken  to  be  zero  and  then  the  printed  stresses  are  additions  to  the 
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initial  state  (i.e.,  superposition  is  valid).  However,  if  the  bounding  surface  model 
is  used,  an  accurate  initiation  of  the  stress  state  is  extremely  important. 


ft  must  be  remembered  that  and  are  negative  when  compressive, 
while  h  is  positive  in  compression.  The  izuts  of  h  are  those  of  stress.  The 
question  of  whether  h  represents  total  or  excess  pore  water  pressure  is  discussed 
in  Section  A3. 

The  initial  states  are  described  by  means  of  simple  linear  equations  in 
depth  y(z),  thus  the  coefficients  a^  ♦  e2  are  dependent  on  the  location  selected 
for  the  origin  of  the  coordinates. 

A  field  example  is  shown  in  Figure  3.  When  a  value  of  g  =  1  is  prescribed 
in  Section  A3,  the  analysis  is  in  term  of  "total"  pore  water  pressure.  The 
vertical  and  horizontal  soil  stresses  can  be  simply  calculated  from  the  soil  weight 
and  a  assumed  constant  value  of  kQ  =  .45.  The  coordinate  system  is  located 
as  shown,  and  psi  units  are  used.  The  following  coefficients  are  found: 
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The  specifications  of  the  initial  void  ratio  e.  and  the  preconsolidation 
pressure  (positive  in  compression)  PQ  are  only  necessary  if  the  bounding  surface 
model  is  used;  they  are  "internal  variables"  for  that  theory  [7]. 

B4.  Node  Point  Array 

The  program  incorporates  two  data  generation  routines  to  assist  the  user 
in  defining  the  locations  of  the  system's  node  points.  The  use  of  these  options 
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can,  lor  example,  enable  one  to  describe  the  locations  of  the  nodes  of  an 
arbitrarily  large  grid  with  as  few  as  five  cards.  Note  that  not  all  numbers 
between  1  and  the  maximum  node  number  NPT  need  correspond  to  actual  nodes 
in  the  body.  For  example,  the  numbering  scheme  shown  in  Figure  4  is  permissible; 
the  coordinates  of  the  non-existent  nodes  15  and  21  need  not  be  specified.  This 
feature  facilitates  the  use  of  the  available  node  point  and  element  generation 
options.  If  the  location  of  a  node  is  prescribed  more  than  once  in  the  input 
and  the  locations  are  not  in  agreement,  the  last  description  is  used.  However, 
if  in  a  second  or  later  description  the  node  number  is  entered  as  negative,  then 
the  previous  location  is  used.  The  utility  of  this  option  is  illustrated  later. 

The  straight  line  or  circular  arc  coordinate  generation  option  may  be  used 
whenever  several  sequential  node  points  lie  along  a  straight  line  or  circular  arc. 
If  such  a  situation  exists,  it  is  necessary  only  to  enter  the  coordinates  of  the 
initial  and  final  points  of  the  sequence  (denoted  by  N*  and  N,  respectively),  and 
the  values  of  INC  and  D.  The  constant  INC  represents  the  difference  between 
any  two  successive  node  numbers  in  the  sequence,  and  D  defines  the  ratio  of 
the  distances  between  any  two  adjacent  pairs  of  points. 

If,  for  a  node  N,  INC  h  0,  intermediate  node  points  are  generated  along 
a  straight  line  (XC  =  YC  =  0)  or  a  circular  arc  (XC  k  0  and/or  YC  £  0)  between 
node  N  and  the  point  described  on  the  preceeding  node  specification  card  N1. 
That  is,  the  coordinates  of  the  points  hP  +  INC,  N*  +  2*INC,  ...»  N  -  INC  are 
each  automatically  found.  For  the  case  of  a  circular  arc  (flagged  by  the 
condition  XC  £  and/or  YC  4  0)  it  is  assumed  to  pass  through  the  end  points  of 
the  sequence  N1  and  N,  and  the  additional  intermediate  point  with  coordinates 
(XC,  YC).  This  intermediate  point  need  not  be  a  node.  The  node  N  for  which 
the  specified  non- zero  value  of  INC  triggers  the  generation  of  the  line  N*  -  N 
can  also  serve  as  the  initial  point  of  a  line  generated  between  it  and  the  point 


described  on  the  next  card.  Thus  the  exterior  nodes  in  Figure  4  can  be  generated 
with  5  cards  for  nodes  1,  6,  42,  37,  and  1,  with  the  last  4  having  specified 
values  for  INC.  If  in  the  second  specification,  node  1  is  entered  as  negative, 
then  its  coordinates  do  not  need  to  be  repeated. 

The  end  points  of  the  sequence  may  be  entered  in  either  order.  For 
example,  the  segments  illustrated  in  Figure  5  could  be  defined  by  specifying  the 
nodes  in  either  the  order  7  -*•  22  (INC  =  5,  D  =  2.0)  or  the  order  22  ♦  7 
(INC  =  -5,  D  =  .5).  The  spacing  of  the  intermediate  points  (nodes  12  and  17) 
is  controlled  by  the  spacing  ratio  D.  A  value  of  D  =  1.0  would  result  in  equally 
spaced  nodes. 

The  interior  node  point  generation  option  automatically  (i.e.  without  any 
prompting  by  the  user)  locates  all  interior  nodes  whose  coordinates  have  not 
been  established  through  the  options  cited  above  (that  is,  all  points  left  undefined 
after  the  input  of  Section  B4  has  been  processed).  The  locations  of  these 
undefined  interior  nodes  are  computed  by  means  of  the  "Laplacian  -  Isoparametric" 
grid  generation  scheme  described  in  [3],  See  also  Section  A2.  Note  that  all 
boundary  nodes  must  be  directly  or  indirectly  specified  by  the  input  to  Section 
B4. 

Figure  6  illustrates  two  grids  that  have  been  prepared  with  the  aid  of 
the  grid  generation  schemes.  Grid  1  was  developed  by  using  the  straight  line 
generation  routine  to  specify  the  exterior  (boundary)  nodes  (only  five  cards  were 
needed  in  Section  B4).  Grid  2  was  developed  in  a  similar  manner,  except  that 
the  straight  line  generation  option  was  also  used  to  define  the  nodes  lying  along 
the  material  interface,  i.e.,  line  3  -►  33.  In  the  generation  of  line  3  *►  33  the 
coordinates  of  nodes  3  and  33  previously  specified  by  line  generation  1  6  and 
31  *►  34  were  used.  This  was  accomplished  by  entering  them  as  -3  and  -33 
(with  blank  x,y  coordinates). 
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B5.  Element  Array 

The  material  number  MN  and  the  initial  state  number  ISNO  must  correspond 
to  the  appropriate  descriptions  given  in  Sections  B2  and  B3.  At  the  time  the 
information  for  the  initial  state  ISNO  is  used  to  initialize  h  for  the  element, 
it  is  also  used  to  initialize  h  for  the  four  nodes  describing  the  element  in 
question.  If  different  initial  states  were  prescribed  for  elements  "a"  and  "c" 
of  Figure  6  and  they  give  initial  values  for  h  at  the  common  nodes  2  and  7, 
which  were  not  in  agreement  the  values  obtained  from  the  element  of  higher 
number  would  prevail.  Because  in  practice  h  is  continuous  such  ambiguous 
situations  should  not  often  arise. 

If  the  body  can  be  divided  into  layers  of  elements,  and  if  the  material 
and  the  initial  state  numbers  MN  and  ISNO  are  the  same  for  several  elements 
within  a  layer  and  for  several  layers,  the  node  numbers  of  these  elements  can 
be  simply  established  by  means  of  the  element  data  generation  option.  To 
generate  a  sequence  of  elements  within  a  single  layer,  node  points  are  specified 
for  the  first  element  only,  together  with  appropriate  values  for  NMIS  and  INC. 

For  example,  the  bottom  row  of  elements  in  the  first  grid  of  Figure  6 
could  be  established  by  entering  either  the  node  numbers  of  element  "a"  and 
the  values  NIMS  =  4  and  INC  =  5  or  the  node  numbers  of  element  "b"  and  the 
values  NMIS  =  4  and  INC  =  -5.  Similarly,  the  left-most  column  of  elements 
could  be  established  by  entering  the  node  numbers  of  element  "a"  and  the  values 
NIMS  =  3  and  INC  =  1,  etc. 

If  several  layers  of  elements  are  of  the  same  material,  and  have  the 
same  initial  state  description  it  becomes  possible  to  carry  this  option  one  step 
further.  For  example,  the  entire  grid  for  the  first  example  of  Figure  6  could 
be  established  by  entering  the  node  numbers  of  element  "a”  and  the  values: 


NMIS 


4 


INC  =5 
NMISP  =  3 
INCP  =  1 

or  alternatively, 

NMIS  =  3 
INC  =  1 
NMISP  =  4 
INCP  =  5 

Grids  generated  with  these  two  specifications  would  differ  only  in  the  order 
that  the  elements  were  numbered. 

Hence,  under  "ideal"  conditions,  the  element  array  for  an  entire  body  can 
be  defined  with  only  a  angle  card  in  Section  B5. 

B6.  Node  Point  Specification  Array 

Boundary  or  interior  node  point  displacement  and  load  specifications  may 
be  given  in  terms  of  either  x-y  (r-z)  components  (8n  =  0  in  Section  B6)  or  local 
Xj-x2  components  (0n  /  0),  see  Figure  7.  If  0  =  0,  the  subscripts  1  and  2  in 

Section  B6  refer  to  x<r)  and  y(z)  (and  thus  IF|  =  IFX,  etc.)  and  if  0n  £  0  they 
represent  Xj  and  x?  (and  thus  IFj  =  1FX  ,  etc.). 

For  each  of  the  two  coordinate  directions,  one  may  specify  the  history 
of  either  a  displacement  (IF  =  1)  or  a  load  (IF  =  0)  by  setting  V  equal  to  the 
magnitude  of  the  applied  quantity  and  IH  equal  to  the  appropriate  history  function 
number  of  Section  Bl.  Specified  displacements  and  loads  are  considered  to  be 
positive  when  they  have  the  same  sense  as  the  positive  coordinate  directions. 
In  addition,  either  the  history  of  the  water  flow  Q  or  the  pore  water  pressure 
h  may  be  specified  by  giving  appropriate  values  for  IHj,  IFj  and  Vj. 
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Figure  7.  Pressurized  Boundary  and  Example 
of  Local  Coordinates. 


V. 


If  several  nodes  in  a  sequence  (N,  N  +  INC,  N  +  2*INC,  .  .  .  NMNC,  N*) 

all  have  the  same  "node  point  specifications",  by  supplying  the  appropriate  values 

for  N,  W  and  INC  they  can  all  be  generated  with  a  single  card.  A  node  point 

may  have  more  than  one  specification  as  long  as  the  accumulative  effect  on 

the  system  of  equations  is  correct.  For  the  same  point,  if  one  specification  is 

a  force  and  the  other  a  displacement  both  in  the  x  direction,  displacement 

specification  will  prevail.  If  both  specifications  were  displacements,  the  result 

would  be  the  second.  Multiple  specifications  which  tend  to  cause  difficulty, 

and  thus  should  be  avoided,  are  those  which  involve  non-zero  values  of  0  . 

n 

Multiple  specifications  at  corners  are  convenient,  and  should  not  be  discouraged. 

Using  the  generation  option  a  uniform  or  linearly  varying  pressure  may 
be  specified  along  a  straight  or  curved  boundary  (or  an  interior  line).  The 
quantities  IFp  IF2,  Vj,  Vj  and  0n  are  left  blank  and  appropriate  values  for  N, 
N*,  PN  and  P^  are  supplied.  For  example,  to  specify  the  boundary  loading 
shown  in  Figure  7,  the  user  would  enter,  on  a  single  input  card,  the  values: 

N  =  11 

=  0. 

n 

N'  =  2 

nc  =  -3 

PN  =  100.0 

PN'  =  50.0 

For  pressure  specifications,  the  points  N  and  N*  must  be  given  in  a 
counter-clockwise  order  if  they  lie  on  an  exterior  boundary  and  in  a  clockwise 
order  if  they  lie  along  an  interior  boundary  (or  "hole").  In  general,  pressure 
specification  cards  should  preceed  all  other  node  point  specifications  in  Section 


If  for  a  given  node,  IF.  =  0,  V.  =  0,  =  0  and  P^,  =  0  it  need  not 

-and,  for  economy,  should  not  —  be  included  in  the  node  point  specification 
array. 


B7.  Solution  History  Segment  Information 

The  analysis  is,  in  general,  time  dependent  due  to  the  consolidation  process 
and  the  history  dependence  of  the  bounding  surface  plasticity  model.  For  a 
non-flow  problem  (IFLOW  =  0)  the  actual  rate  at  which  time  passes  is  not 
important  (because  the  bounding  surface  model  is  rate  independent),  however, 
for  the  purpose  of  modeling  the  history  effects  it  is  still  convenient  to  think 
in  terms  of  time.  For  a  linear  elastic,  non-flow  problem  the  only  role  of  time 
is  to  represent  the  loading  history;  if  only  final  results  are  desired  then  only 
one  time  step  is  required. 

For  convenience,  the  solution  history  is  broken  into  one  or  more  history 
segments.  One  card  is  required  in  Section  B7  for  each  segment.  The  time  at 
the  end  of  a  given  segment  is  denoted  as  TIME;  it  is  assumed  that  the  first 
segment  begins  at  t  =  0.  The  number  of  time  steps  into  which  a  given  segment 
is  to  be  divided  is  prescribed  as  NMIS.  Within  a  given  ’history  segment”  the 
ratio  of  two  successive  time  steps  is  equal  to  the  prescribed  spacing  ratio  D; 
a  value  of  1.0  gives  equal  time  steps.  The  role  of  the  time  step  spacing  ratio 
is  analogous  to  the  length  spacing  ratio  used  in  Section  B4  and  Figure  5. 

The  selection  of  appropriate  time  step  lengths  is  complicated  by  the  fact 
that  two  distinct  processes  are  involved,  i.e.  water  flow  and  soil  plasticity. 
Thus  a  certain  amount  of  experimentation  with  successively  smaller  time  steps 
will  often  be  required.  In  this  process  several  factors  should  be  considered. 
Abrupt  changes  in  step  size  should  be  avoided  (judicious  use  of  the  spacing  ratio 
D  can  facilitate  smooth  transitions  from  small  to  large  time  steps,  etc).  Abrupt 
changes  in  applied  loads  or  displacements  will  cause  large  flow  gradients  and 
require  small  time  steps.  Further  information  concerning  step  size  for  flow 
problems  is  to  be  found  in  [1,  21.  A  certain  amount  of  oscilation  of  the  solution 
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is  to  be  expected  and  usually  can  be  tolerated.  A  minimum  of  10-20  steps  are 
usually  required  by  the  bounding  surface  plasticity  model  in  proceeding  from  a 
nearly  hydrostatic  stress  state  to  failure  conditions. 

C.  End  Card 

The  function  of  this  card  is  to  signal  the  end  of  the  problem;  the  program 
then  proceeds  to  the  next  stacked  job  (if  any). 
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V. 


EXAMPLE 


The  following  example  is  intended  to  illustrate  the  program  input  and 
output;  the  accuracy  of  the  results  is  discussed  in  [1  ].  The  problem  configuration 
is  shown  in  Figure  8.  It  is  an  idealized  representation  of  the  consolidation  of 
a  soil  mass  that  is  free  to  drain  both  from  the  top  and  a  central  sand  drain. 
The  input  file  and  the  output  for  the  first  3  time  steps  are  given  in  the  Appendix. 
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4 

392E*00 

16 

3 

303E-02 

-3.  238E-02  -1. 

733E-02 

1 

043E02  1 

101E+01 

-6 

476E+00 

17 

-7 

983E  03 

3  151E-04  -7 

130E-04 

-4 

022*  03  -1 

397E*00 

1 

030E-01 

18 

-1 

307 t -02 

-7.  8T1E-03  -2 

171E-03 

-2 

633E-03  -2 

613E+00 

-1 

568E+00 

19 

-l. 

100t-02 

6.  083E-03  -2 

883E-03 

2 

910E-03  -2 

200E*00 

i. 

217E+00 

20 

-1 

244E-02 

-1. 374E-02  -2 

688E-03 

3 

S76E  03  -2 

439E+00 

-2. 

749E+00 

21 

-1 

237b-02 

1.400E-02  -2 

294E-03 

4 

708E  03  -2 

314E+00 

2. 

799E+00 

22 

-1 

320E-02 

-2.  163E-02  -1 

812E-03 

6 

197b  03  -2 

640E-00 

-4 

327E+00 

23 

-1 

677E-02 

2.  337E-02  -1 

673E-03 

3 

432E-03  -3 

333E*00 

1 14E+00 

24 

-1 

653E-02 

-3.  019E-02  -2 

707E-03 

-2 

890E-03  -3 

306E+00 

-A 

038E*00 

23 

7. 

290t  03 

8  349E-04  -6 

083E-04 

-3 

48  IE -03  1 

438E+00 

1. 

670E-01 

26 

1. 

394E-02 

-6  324E-03  -1 

426E-03 

-3 

803E  04  2 

788E+00 

-1 

263E-00 

27 

1. 

341E-02 

6.  662E-03  -1 

717E-03 

9 

303E-04  2 

691E+00 

1 

332E-00 

28 

1 

364E-02 

-1  339E-02  -1 

749E-03 

2 

386E-03  2 

729E-00 

-2 

719E*00 

29 

1 

286E-02 

1  31 BE -02  -1 

397E-03 

3 

470E  03  2 

372E*00 

3 

036E+00 

30 

1 

194E-02 

-2  1 IOC— 02  -1 

474E  -03 

3. 

243E  03  2 

389E-00 

-4 

219E-00 

31 

1 

330E-02 

2.  331E-0?  -1 

692E-03 

8. 

303E  04  2 

639E+00 

3 

I02E-00 

32 

1 

424E-02 

-3  231E-02  -2 

261E  03 

2 

113E-03  2 

B47E+00 

-6 

303E-00 

33 

-3 

392E  03 

1  410E-03  -3 

9B9E-03 

-3 

608E  04  -6 

783E-01 

2 

821E-01 

34 

-3 

879E-03 

-6.  070E-03  -2 

131E-04 

-8 

804E-04  -1 

176E*00 

-1 

214E*00 

33 

-4 

919E  03 

7.  397E-03  -3 

924E-04 

6 

068E  04  -9 

838E-01 

1 

303E-00 

36 

-3 

263E-03 

-1  326E-02  -4 

291 E— 04 

372*  03  -1 

033E*00 

-2 

632E*00 

37 

-3 

402* -03 

1.  360E-02  -4 

130E-04 

7B6E-03  -1 

OBOE +00 

3 

120E*00 

38 

-3 

327t  03 

-2  113E-02  -4 

337E-04 

1 

497t  03  -1 

103E-00 

-4 

230E*00 

39 

-6. 

236E-03 

2  306E-02  -3 

336E-04 

190E  03-1 

231E»00 

3 

012E-00 

40 

-3 

322E-03 

-3  168E-02  -7 

683E-04 

-1 

174E-03  -1 

064E*00 

-6 

336E*00 

41 

2 

380E-03 

1  412E-03  -1 

064E-04 

-9 

671E  04  3 

160E-01 

2 

B23E-01 

42 

3 

013E  03 

-3  393C-03  -2 

33IE-04 

-8 

283E  03  1 

003E+00 

-1 

079E*00 

43 

4. 

902E -03 

7  763E-03  -3 

223E-04 

2 

484E  04  9. 

eosE-oi 

1. 

3S3E*00 

44 

3 

107E  03 

-1.277E-02  -3 

633E-04 

6 

786E-04  1 

021E*00 

-2 

333E«00 

43 

4. 

923E-03 

1  337E-02  -3 

813E-04 

1 

026E-03  9 

845E-01 

3 

U3E*00 

46 

4 

916E-03 

-2  1 12E-02  -4 

104E-04 

7. 

384E-04  9 

831E-0I 

-4 

223E+00 

47 

3 

627b-03 

2  499E-02  -4 

938E-04 

2 

740E-04  1 

125E+00 

4 

997E+00 

48 

3 

307t-C3 

-3.  224E-02  -6 

USE -04 

1 

131E  03  1 

101E+00 

-6 

447E«0O 

49 

-1. 

143E-03 

1 .  394C-03  2 

030E-03 

2 

418E  04  -2 

290E-01 

3 

I88E-01 

30 

-2 

031E-03 

-3.3W4E-03  l 

326E-03 

-2. 

347E  04  -4 

062E-01 

-l 

065E+00 

51 

-1 

676E-03 

8  120E-03  -2 

471E-03 

1 

018E-04  -3 

331E-01 

1 

624E-00 

32 

-1 

673E  03 

-1.27CE-02  -4 

311E-05 

3 

84 IE  04  -3 

346E-01 

-2. 

340E-00 

33 

-1 

343E-03 

l  373E-02  -6 

298E-0? 

3 

1I4E  04  -3 

090E-01 

3 

146E»00 

34 

- 1 

498E  03 

-2.  120E-02  -8 

441E-0S 

4 

477b  04  -2 

997E-01 

-4 

240E*00 

33 

-1 

636E-03 

2  492F-02  -1 

123E-04 

3. 

043E  04  -3 

31  IE-01 

4 

983E*00 

36 

-9 

618E-04 

-3.  205E-02  -l 

681E-04 

-7 

464E  04  -1 

924E-01 

-6 

409E*00 

37 

4 

406E  04 

1.  693E-03  -2 

937E-05 

-4 

079E  04  8 

81  IE-02 

3 

209E-01 

38 

9 

162F-04 

-3  19CE-03  -6 

103E-03 

4 

442*  03  1 

832E-01 

-1 

038E+00 

39 

9. 

984E-04 

8  224E-03  -6 

636E-03 

3. 

220E  OS  I 

997E-0I 

1 

643E*00 

60 

1 

130E  03 

-1  263E-02  -7 

332E-OI 

2 

082E -  06  2 

260E-01 

-2 

329E-00 

61 

1 

182F  03 

1.  382E-02  -7 

87BE-03 

2 

333E  04  2 

364E-01 

3 

145E*00 

62 

1 

298E-03 

-2  134E-02  -8 

63JE-03 

-7 

528E  03  2. 

S96E-0I 

t4 

268E+00 

63 

1 

33SE  03 

2  307E-02  -1 

039E-04 

2. 

846E  06  3 

116E-01 

015E*00 

64 

374E-03 

-3. 23BE-02  -1 

049E-04 

8 

426E  O'.  3 

147E-01 

-6 

476E+00 

TAO-RZ 
I  991E+00 
1  088E+00 
I  4S9E-01 
7  407E-02 
1  493E-0I 
4  2B3E-0I 


-6 

-7 

-I 

-3 

-3 

-3 

-3 

-3 

-2 

-3 


S40E+00 
631E+00 
732E+00 
322E+O0 
492E+00 
390E+00 
29 1 E  +  DO 
097E*00 
781E+00 
471E+00 
1  426E-01 
4  343E-01 
3  770E-01 

3  376E-01 

4  3B7E-01 
3  623E-01 
3  349E-01 

415E-01 
2I7E-01 
832E-01 
434E-01 
497E-01 
193E-01 
948E-0I 
384E-0I 
323E-01 
979E-03 
263E-02 
848E-02 
383E-02 
261E-02 
674E-02 
1  067E-01 

1  337E-01 

2  I29E-02 
061E-02 
431E-02 
310E-02 
630E-02 
208E-02 
876E-02 
224E-01 
061E-03 
031E-03 
94  IE-03 
021E-03 
260E-02 
688E-02 
246E-02 
362E-02 
874E-03 
222E-02 
33 IE-02 
306E-02 
576E-02 
731E-02 
077E-02 
098E-02 


45 


7  1B8E*01 
7.  212E+01 
7  392E+01 
7  23 IE* 01 
7  769E-01 
7  079E*01 


093E-01 
397E*00 
368E*00 
739E-0I 
036E-01 
313E-0I 
239E-0I 
983E-01 
l  373E»00 

1  04 3E+00 
-4  022E-01 
-2  633E-01 

2  910E-01 

3  376E-01 
708E-0I 
197E-0I 
432E-01 
B9OE-0I 
481E-0I 
803E-02 
305E-02 
386E-01 
470E-01 
243E-0I 

8  303E-02 
2  II3E-01 
-3.  608E-02 
-8  804E-02 
6  068E-02 
l  372E-0I 
1  786E-01 
I  497E-0I 

1  190E-0I 
■1.  174E-01 
-9  671E-02 

283E-03 
484E-02 
786E-02 
02&E-01 
384E-02 
740E-02 
t.  151E-0J 

2  418E-02 
-2  347E-02 

1  018E-02 
3.  84 IE-02 

3  114E-02 

4  477E-02 

3  043E-02 
-7  464E-02 
-4.  079E-02 

4  442E-03 
3  220E-03 

2  0B2E-04 
2  333E-02 

-7.  328E-03 
2  846E-04 
8  426E-02 


184E«01 
968E-01 
072E+02 
043E-02 
082E-02 
036E+02 
I03E*02 
02 IE 402 
133E+02 
02 IE 402 
679E+01 
386E+01 
S03E«0I 
443E+01 
93SE+01 
33BE*01 
I  0I3E+02 

8  978E+0I 
026E*02 
CI1E402 
038E+02 
934E40I 
036E+02 
760E*0l 
077E»02 
3B3E-01 
921E»0I 

9  775E+01 
1  004E+02 
9.  629E+01 
1.  018E-KJ2 

473E«01 
034E+02 
231E+01 
012E+02 
992E*0I 
023E-02 
9.  83GE+01 
1  041E+02 
9  648E+0I 
I.  061E*02 
9  433E+01 
000E-02 
B39E+01 
0l3E*02 
7I2E+01 
027E-02 
343E+0I 
044E+02 
321E*01 
005E+02 
9I8E40I 
0I9E+02 
767E*01 
034E*02 
386E*0I 
032E+02 
364E*0l 


U-l 


j-'-A 


NODE  DlSPt-ACFHENlS 
NO  V  V 

1  -1  63SE-19  -1  4266-19 

2  -3  78CE-19  -I  43OE-02 

3  -4  026E-19  -1  29OE-02 

4  -3  994E-19  -fe  463E-03 

3  -4  021F-19  -»  B39E-02 

4  -3  966E-19  -7  612C-03 

7  -4  Q11F-19  -3  641E-02 

8  -3  911C-19  -I  99IE-07 

9  -2  069E-19  -I  2B6E-01 
10  -1  429E-19  -9  003E-19 
U  -3  696E-02  -4  204E-03 

12  -3  09 IE  02  -1  377E-02 

13  -3  042F-02  -4  172E-03 

14  -3  029E-02  -l  987E-02 
13  -4  933E-02  -4  780E-03 

16  -4  7*04-02  -3  333E-0P 

17  -4  433E-C2  -I  677E-02 

18  -6  278E-02  -3  387E-02 

19  3  078E-20  -2  012E-18 

20  4  419E -03  9  723E-04 

21  -2  210E-03  -8  856E-03 

22  -1  213E  03  -2  464E-03 

23  2  166E-04  -1  632E-02 

24  8  838E  04  -3  022E-03 

23  3  233E-03  -2  353E-02 

26  3  140E-03  -1  694C-04 

27  -2  147t-03  -2  383E-02 

28  -7  323E-20  -2  99SE-1S 

29  -1  133E-C2  3  771E-03 

30  -1  237E-02  -3  797E-03 

31  -l  303E-02  -3  2B7E-03 

32  -1  283E-02  -1  346E-02 

33-1  120E-02  1  04CE-03 

34  -1  106E-02  -1  972E-02 

33  -1  408E-02  6  039E-03 

36  -1  398E-02  -2  863E-0= 

37  4  046E-20  -4  016E-18 

38  3  033E-03  1  612E-03 

39  9  277E-04  -3  181E-03 

40  4  S97t-04  2  369E-03 

41  9  432F-04  -1  138C-02 

42  7  417E-04  4  478E-03 

43  8  821E  04  -1  493C-02 

44  3  720E-04  8  291E-03 

43  -2  164E-03  -2  202E-02 

46  -3  790F-20  -3  00 IE- 18 

47  -3  731E-03  l  209E-03 

48  -4  046E-03  -4.  138E-03 

49  -4  404E-03  3  366E-03 

50  -4  723E-03  -9  40OE-03 

31  -4  396E-Q3  3  943E-03 

52  -i  034E-03  -1  493E-02 

33  —  j  025E-03  9  943E-03 

34  -o  219E-03  -2  310E-02 

35  2  131E-20  -6  012E-18 

36  1  409E-03  1  616E-03 

37  P  20 IF -04  -3  827E-03 

38  3  343E-04  4  196E-03 

57  J  323F.-04  -8.  372E-03 

60  :  733E-04  7  218E-03 

61  2  276E-04  -l  414E-02 

62  -3  I81E-C3  1  096E-02 

63  -1  192E-03  -2  047E-02 

64  -1  19JF-20  -7  00CE-I9 
63  -8  811E  04  l  372E-  03 

66  -9  313E-04  -3  633E-03 

67  -l  046E-03  4  383E-03 

68  -l  214E-03  -8  043E-03 

69  -1  130c -03  7  626E-03 

70  -1  44«,e-C3  -1  34IE-07 

71  -I  6t>9£-03  1  132F-02 

72  -1  478E-03  -2  133E-02 

73  4  0l7fc-lB  -3  839E-18 

74  7  97/E— 18  1  637E-03 

73  8  027t-18  -3  339E-03 

76  7  963E-18  4  690E-03 

77  8  024E-18  -7  973E-03 

78  7  930c- 18  8  00 1C- 03 

79  8  022E-1B  -1  364E-02 

80  7  933E- 18  1  178E-02 

■1  3  731C-18  -2  03 1 E- 02 


■  TIME  0  000 1 J"  NO  ITERATION  WAS  REQUIRED 


PMENT  El  EMENT  STRAINS  AND  STRESSES 

NO  EPSILON-R  EPSILON-Z  E’S-TNETA  OAMMA-RZ  SIONA-R  SIOMA-Z  SIC-THETA  TAU-RZ  h 


1  -7  471E-02  -1  7236-02  -3  202E-02  -3  294E  02  -1 

2  -1  430E-01  -3  8996-03  -A  213E-02  1  270E02  -2 

3  -I.  405E-0I  1.  1A4E-02  -A  023E-02  1  789E  03  -2 

4  -1.411E  01  -2  274E-02  -6.  047E-02  -I  3S4E  Ol  -2 

3  -1  3»«  01  2.  2846-02  -9  9796-02  4  816E  03  -2 

A  -1  376E-01  -9  1A3E-02  -9  909E-02  3  2B4E-04  -2 

7  -1  2686-01  3  400E-02  -9  4336-02  1  979E  02  -2. 

8  -1  1846-01  -1  2A8E-01  -9  083E-02  8.  2A3E  02  -2 

9  4.  4946-02  -3  785E-03  -1  S04E  02  -3.WKE-02  8 

10  8  2S8E-0?  -1  3AAE-02  -3  09BE-02  7  989E  03  1. 

11  7622602  1.  479E-02  -3  081E-02  A  279C-03  1. 

12  7  817E  -02  -2  878E-02  -3  038E-02  3  989E-03  1. 

13  7. 874E-02  3  2046-02  -2  9926-02  68236-03  1 

14  8  0896-02  -3  94AC-02  -2  817E-02  9714E-03  1. 

19  7.  740E  02  A  1 92C-02  -2  493E-02  2  799E-02  1 

1A  A  794E-02  -9.  3126-02  -2  493E-02  2.  81AE  02  1. 

17  -A  9966-03  1  972E-03  -1  30AE-03  -8.  249E-03  -1. 

18  -1.2036  02  -1  SI  SC-02  -4  2966-03  -4  3896-03  -2 

19  -1  14AE-02  1  389E-02  -9.  93AE-03  4  962603  -2. 

20  -1.  319E-02  -3  000E-02  -9.  233E-03  3  921F  03  -2. 

21  -1  338E-02  3  923E-02  -4.  A37E-03  7  0A4E-03  -2. 

22  -1.  493S-02  -3  AA2C-02  -3.  720E-03  1.  131E-02  -2. 

23  -1  839E  02  7  282P-02  -3  132E-03  1  0O4E-02  -3 

24  -1.  3936-02  -9.  A42C-Q2  -4  A42E-03  -2  737t-03  -3. 

29  7.  4A1E-03  1 .  7AOC-03  -1  003E-03  -A  1329  03  1 

2A  1  A30t-02  -1  230E-O2  -2.  497E-03  -1  421E  03  3 

27  1  7A7t-02  1  419E-02  -3.  OAAE-03  1.  389C-03  3. 

28  1.  740E-02  -2.  930E-O2  -3.  137E-03  4.  407t-O3  3. 

29  1  A33E-02  3  A43E-02  -2.  8A2E-03  9.  972P-03  3. 

30  1  939E-02  -3  429C-02  -2  394E-03  9  8O9E-03  3 

31  1.  903E  03  7.  19 SC- 02  -2  731E-03  3.  494E-03  3 

32  1.  487E  02  -9  839E-02  -3  470E-03  3  498E  03  2 

33  -1  946E  03  2  702C-03  -l.  A92E-04  -1  48AE  03  -3 

34  -4  0426-03  -1  171E-02  -9  48SE-04  -1  130E  03  -8. 

33  -3  861 E- 03  l  3A7E-02  -8.  307E-04  8  923E-04  -7. 

3A  -3  8246-03  -2  8636-02  -9  320E-04  2  AllE  03  -7 

37  -4  41 AE  03  3.  703E-O2  -8  8026-04  3  4436-03  -8 

38  -4.  910E-03  -3  417E-02  -8  323E  04  3  097t  03  -9 

39  -4.  8196-03  7  1A0E-02  -9  91AE-04  2  038E  03  -9 

40  -3  080E-03  -9  82 IE- 02  -1  3886-03  -1  032E03  -A 

41  1  8816  03  2  829E-03  -1  444E-04  -1  49BE  03  3 

42  4.  A3AE-03  -1  080E-02  -3  932E-04  -7  4A9E  01  9 

43  9  499E-03  1  A34E-02  -9  4716-04  4B81E01  1 

44  3  483E-03  -2  B33E-02  -A  U7E-04  l  33AE03  1 

49  3  393E-03  3.  728E-02  -A  309E-04  1  738E-03  1 

4A  9  3466-03  -3  412E-02  -A  9776-04  1  A9AE  03  l 

47  3  479E-03  7.  13AE-02  -7  922E-04  1010603  1. 

48  3  2406-03  -9.  697E-02  -9.  39AE-04  1  288E  03  1 . 

49  -4.  A13E  04  3  0A1E-03  -1.301E  09  -2  819E-09  -9 

90  -1  0506  03  -1  092E-02  -9.  929E-03  -3.  808E-09  -2. 

91  -9  9A3E  04  1 .  A84E-02  -1.136E-04  1  012E-04  -1 

92  -7. 379E-04  -2  7/AE-02  -1  941E-04  A  3286-01  -1 

33  -7  1196  04  3  747E-02  -1  7346-04  8  4A2P  01  -I 

34  -7  1 13E  04  -3  417E-02  -2  0016-04  7  OllE  01  -1 

93  -4  807t  04  7  125E-02  -2  321E-04  3  0946  01  -9 

3A  A  981E-04  -9  87AE-02  -3  383E-04  -t,  089E-09  l 

97  3  1 SPE-04  3  09CE-03  -2  101E-03  -4  27CE  04  A 

98  8  B40E-04  -1  027E-02  -9  B93E-03  -2  8aOE  01  l 

99  I  21AE-03  1  A94E-02  -8  106E-09  1  04AE  04  2 

AO  1  sso^-oa  -2  733E-02  -9  2036-03  1  217t  04  2 

61  1  483S  03  3  791E-02  -9  886E-03  2  921E  01  2 

62  1  AS AC -03  -3  431E-02  -1  1046-04  A  039E  03  3 

«3  1  879E-03  7  147C-02  -1  2936-04  2  4  186  0"  3 

A4  1  830E-03  -9  9406-02  -1  230E-04  A.  339E  04  3 


4946+01  -3  447E+00  -A  4046+00  -3  2546*00  3.  S13E+01 

9006*01  -7  798E-01  -l  2436*01  1  2706+00  9  3226*01 

8UE+01  2.  327E+00  -1.  2096*01  1.  7B9E-01  5  932E+01 

622E+01  -4  349E+O0  -1.209E+01  -1  394E-02  9.  404E+01 

788E+01  4  3A7E+00  -1  1936*01  4.  81AE-01  A.  238E+01 

7966+01  -1  033E+01  -1 1816*01  3  2846-02  9  077E*0l 

33AE+01  A  801E+00  -1  0876+01  1  9796*00  A.  83SE+01 

3726*01  -2  93AE+01  -1  0176*01  8  2A3E+00  9  4426+01 

9096+00  -7  370E-01  -3  0076*00  -3  9046*00  1  087E+02 

A1AE*01  -3  1326+00  -A  1196+00  7989E-01  1  0346*02 

3246*01  2  998E+00  -A  1A3E+00  A  2736-01  1  124E+02 

S63E+01  -9  737E+00  -A.  07AE+00  3  989E-01  1  0376*02 

9736+01  A  4086*00  -5  9046+00  A  8296-01  1  1736*02 

A17E*01  -1.  189E+01  -3  A34E  +00  9  714E-01  9  892E«01 

318E+01  1  23OE*01  -4  983E*00  2  799E+00  1  2A4E*02 

331E+01  -1  908E+01  -4  987E+00  2  81AE+00  9.  071E*01 

393E+C0  3  143E-01  -3  013E-01  -B  249E-01  9  742E«0l 

409E400  -3  090E+00  -8  992F-01  -4  9B9E-01  9  394E+01 

291E-00  2  778E*00  -1  107E*00  4  982E-01  9891E*01 

639E*00  -6  OOOEiOO  -l  047E+00  5  721E-01  9  044E401 

A73E+00  7.  049E*00  -9  274E-01  7  0A4E-01  1  033E«02 

989E-00  -1.  132E+01  -7  441E-01  1  131E+00  8  3A8E-01 

A77E«00  1  49a£-»0I  -A  304E-01  1  004E*00  1  103C*02 

190E+00  -1  928E*01  -9  284E-01  -2.  737E-01  7.  A9OE-01 

492E«00  3.  320E-01  -2  011E-01  -A  132E-01  1.024E«02 

2A1E+00  -2  4A0E+00  -4  913E-01  -1  421E-01  1.  0O4E»02 

334E4O0  2  B9BE-00  -6  133E-01  1.  383E-01  10A0E»02 

«79E*00  -9  900E+00  -6  274E-01  4  407E-01  9.  A83E-01 

3Q9E+00  7  2BAE*00  -3.  724E-01  9.  972E-01  1.103E+02 

079E-»00  -1.  089E*01  -3.  187E-01  3.  909E-01  9.  1A3E»01 

003E+00  1  439E+01  -3  4A3E-01  3.  494E-01  1.170E+02 

979E*00  -1  978E+01  -A  940E-01  3  498E-01  8.  243E*01 

893E-01  9  404E-01  -3  383E-02  -1  48AE-01  1  O02E+O2 

084E-01  -2  342E*00  -l  098E-01  -1.13OE-01  9  A83E-KM 

722E-01  3  134E*00  -1  701E-01  8  923E-02  1022E*02 

A48E-01  -3  726E+00  -1  8A4E-01  2.  A11E-01  9  343E*01 

B32E-01  7.  40AE+00  -1  7A0E-01  3  443E-01  10A2E-02 

B20E-01  -1  08JE+01  -1  709E-01  3  097E-01  8808E*01 

A38E-01  1  432E+01  -I  983E-01  2  038E-01  1  129E*02 

1996-01  -1  9A4E+01  -2  777E-01  -1.  032E-01  7  922E*01 

7A1E-01  9  A98E-01  -2  888E-02  -1  498E-01  1.0UE«02 

313E-01  -2  160E*00  -7.  BA3E-02  -7  469E-02  9.8876*01 

1006*00  3  2696*00  -1  094E-01  4  881E-02  1  0446*02 

097E+00  -9  A07E+00  -1  223E-01  1.  33AE-01  9  934E*01 

0806+00  7.  497E+00  -1  2A2E-01  1  738E-01  1  0846*02 

049E+00  -1  082E+01  -1  319E-01  1  6966-01  89996*01 

09AE+00  1  427E+01  -1  904E-01  1  010E-01  1  1316*02 

01BE+00  -1  979E+01  -1  879E-01  1  28BE-01  8  0776*01 

2^AE-02  A  121E-01  -2  602E-03  -2  819E-03  1  0076*02 

100E-01  -2  104E*00  -1  103E-02  -3  808E-03  9  779E*01 

913E-01  3  3A7E+00  -2  271E-02  1  012E-02  1  032E*02 

919E-01  -3  991E*00  -3  082E-02  A  32BE-02  9429E*01 

424E-01  7  499E*00  -3  4A8E-02  8  462E-02  1  072E*02 

423E-01  -l  083E*01  -4  001E-02  7  011E-02  8  89IE*01 

A13E-02  1  4256*01  -3041E-02  3034E-02  1  138E*02 

39AE-01  -1  979E+01  -A  7A5E-02  -4  0896-02  8  002E*01 

3046-02  A  180E-01  -4  203E-03  -4  270E-02  1  008E*02 

7A8E-01  -2  034E+00  -1  199E-02  -2  8A0E-02  9  822E*01 

432E-01  3  3876*00  -I  A22E-02  1  04AE-02  1  037E*02 

761E-0I  -9  3066*00  -I  841E-02  1  2176-02  9  4796*01 

9666-01  7  3036*00  -1  9776-02  2  9216-02  1  077t*02 

3126-01  -1  0866*01  -2  2086-02  A  0336-03  8  9306*01 

7906-01  1  4296*01  -2  5056-02  2  4186-02  I  1456*02 

A99E-01  -1  9886*01  -2  46AE-02  A  333E-02  8  0186*01 


AT  TIME  0  OOCJf  NO  ITERATION  WAS  REQUIRED 


El  EKENT  Ei  6MENT  STRAINS  AND  STRESSES 

NO  EPSILON-R  EPSILON- Z  EPS-THE I A  CAMMA-RZ  SIOMA-R  SIQMA-Z  SIC-THETA  TAU-RZ  h 

1  -6  A82F-02  -2  3116-02  -2  0646-02  -3  0146-02  - 1 . 3366+01  -*  A81E+00  -3  727E+00  -3  014E+0O  3  771E+01 


2  -I  2966-01  3.  A34E-03  -5  3336-02  1  347t  02  -2 

3  -1  23IP-01  2  2126-03  -3  3A0E-02  2  402E  03  -2 

4  -1  247t-01  -I.  1O4E-02  -3  3466-02  7  3306  04  -2 

3  -1  2136  01  1  028E-02  -3  2O8E-02  A,  82 IE  03  -2 

A  -1  222E-01  -4  1BAE-02  -3  237E-02  -4  A12E  03  -2 

7  -1  173E  01  3  24AE-02  -3  OME  02  1.7346  02  -2 

8  -1  1466-01  -1  JTOE-Ol  -4.  913E-02  7  S38E  02  -2 

9  3  30  76  02  -9  922E-03  -1  3716-02  -3  9606-02  A 

10  A  092F  02  -3  O17E-03  -3  1326-02  1  139E  02  1 

11  3  383E  02  -4  910E-04  -3  I41E-02  A  310E-03  1 

12  3  710E-02  -1  048E-02  -3  08AE-02  3  331E  03  1 

13  3  A81E  02  1  322P-02  -2  9A7E-02  7.  9A3E  03  1 

14  A  034E-0?  -4  A94E-02  -2  877E-02  7  6726  03  1 

13  A  191E  02  A  114E-02  -2  A30E-02  2  9B2E-02  l 

1A  3  370fc-02  -1  2376-01  -2.  728E-02  3  183602  1. 

17  3  438E-03  -2  417E-03  -2  U9E-03  -1.060E  02  A 

18  A  292-— 03  -3  411E-03  -3  4A9E-03  -3  2436  03  1 

19  3  106E  03  —3.  394E-03  -A.  A3AE-03  3  A71F  03  1 

20  4  018S-03  -0  81AE-03  -A  294E-03  A  9A1E  03  8 

21  3  849E  03  1. 33AE-02  -3  6726-03  7321603  7 

22  2  481 E— 03  -4  0/0E-02  -4  A98E-03  1.312E  02  4. 

23  -3  2196-03  7  179E-02  -4  0426-03  1  133E  02  -6 

24  -3  9086-03  -1  2346-01  -A  010E-03  -3  A77fc-03  -7 

23  9  1786  04  -2  0366-03  -8.  884E-04  -3  849E  -03  1 

2A  4  A27t  03  -1  783E-03  -2.  34AE-03  -2  3496-03  9 

27  6  804E-03  -3  010E-03  -3  033E-03  2  1326  03  1. 

28  3  832F  03  -8.  328E-03  -3  089E-03  4  936603  1 

29  3  002F-03  1  409E-02  -2  787E-03  A  718E-03  1 

30  3  443E-03  -3  7226-02  -2  309E-03  A  402F  03  A 

31  4070603  A  978E-02  -2  7AAE-03  3.1636-03  8. 

32  7  7976  03  -1  273E-01  -3  737E-03  2  A39E03  1 

33  2  A1AE-03  -1  247E-03  -2  9B4E-04  -2  24AE  03  3 

34  4.  A47E-03  -7  9376-04  -7  74AE-04  -A  143E  04  9 

33  4  A22E-03  -1  6686-03  -1  10SE-03  8  2336-01  9. 

3A  4  924E-03  -7  3086-03  -1  207E  03  2  949E  03  9 

37  4  136E  03  1  4976-02  -1  130E-03  3  8936  03  8 

38  3  839E-03  -3  731E-02  -1  140E-03  3  199E  03  7 

39  3  449E  03  A  972E-02  -l  3IAE-03  2  292*  03  A 

40  2  793E  03  -1  2A9E-01  -l  730E-03  -2  2486-01  3 

41  -8  380E-04  -8  3846-01  -8  433E-03  -9.  789E  04-1. 

42  -7  0176-01  -7  2636-03  -2  913E-04  -1  014E  03  -l 

43  2  1A0E  04  -8  4A1E-04  -4.  A44E-04  A  394E  04  4. 

44  9  131E-05  -A  8136-03  -3  3186  04  l  493E  03  1 

43  3  A94E-03  1 .  324E-02  -3  397E-04  1871603  7. 

4A  -1  442E-04  -3  731E-02  -3.  934E-04  1  927E  03  -2 

47  4. 409E  O'  A  9396-02  -7  2286-04  7.3116-01  8. 

48  1  973E  03  -1.2736-01  -9  817E-04  4  238E  04  3 

49  1  031 E- 03  -7  4106.-04  -5  8O7E-03  -3  8A3E-04  2 

50  2  004E  03  4  3126-04  -1  4A3E-04  7  097fe  OA  4 

31  2  133E-03  -3.  A49E-04  -2  10AE-04  7.  328E-03  4. 

32  2  408E  03  -A  489E-03  -2  37/6-04  7  032*  04  4 

33  2  4A2P-03  1  33AE-02  -2  814E-04  9  1A4E  04  4 

34  2  344E-03  -3.  747E-02  -3  19JE-04  6323E01  3 

33  2  AA2F  03  A  918E-02  -3  729E-04  8  299E  04  3 

3A  3  0736-03  -1  2736-01  -4  42JE  04  3  SA7fe  04  A 

37  -1  381 *-04  -6  7346-04  9  203E-0A  3  7066  03  -2 

38  -4  947fc-03  7  1106-04  3  298E  OA  -2  140E  04  -9 

39  2  9186-04  -3  4826-04  -1  943E-0S  A  307t  03  3 

AO  4  7139-04  -A.  0376-03  -3  143E-03  1  977E  04  9 

At  3  97/604  1  3126-02  -3  984E  03  2  470E  04  1 

42  8  032*-04  -3  7336-02  -3  334E-03  1  293E  04  1 

A3  1.0926-03  A  9216-02  -7  2826-03  1  3036  04  2 

A4  1  3386-03  -1.2796-01  -8  9236-03  -2  1726  04  2 


5916+01  7  7A7E-01  -1.  1116*01  1.347E+00  A  03SE*01 
S01E*01  4  423E-01  -1.072**01  2  4026-01  A.  2A 1  E*0 1 
49SE*01  -2  208E*00  -1.  0A9E*01  7 .3306-02  A  071E+01 
4306*01  2  03AE+00  -1  042E+01  A  821E-01  A  388E*01 
4446*01  -8.  373E*00  -1.0476*01  -4  A12E-01  3.  731E*01 
347t*01  A  491E+00  -1.  00AE+01  17346+00  71346+01 
293E*01  -3  193E+01  -9.  823E*00  7  3386*00  4  973E*01 
A13E+00  -1  9B4E+00  -3  141E*00  -3  9A0E*00  1  028E*02 
21BE+01  - 1  0036*00  -A.  304E*00  1  139E*00  1  033E*02 
1 1 7E*01  -9  8216-02  -A.  2826+00  6  310E-01  1  047E+02 
142E+01  -2  097E+00  -A  172E+00  3.  331E-01  1  0326*02 
13e£+01  2  A43E+00  -3  9346*00  7  9A3E-01  1  0916*02 
207E+01  -9  3886*00  -3  754E+00  7  A72E-01  9  712E+01 
23BE  +01  1  223E+01  -3.  2A0E+00  2  982E+00  1  2326*02 
114E+01  -2.  373E+01  -3  43AE+O0  3  1836*00  8  209E+01 
91AE-01  -4  834E-01  -4  238E-01  -1  OAOE+OO  1  0036+02 
248E+00  -1  O82E+0O  -1  094E+00  -3  2436-01  9  904E+01 
021E+00  -7  1A8E-01  -1.  331E+00  3  A71E-01  9  83AE+01 
037E-01  -1.  7A3E+0O  -1.  239E+00  A.  9A1E-01  9  801E+01 
A98E-01  2  A13E+00  —1.  134E+00  7.  321E-01  1.0216*02 
9A2E-01  -8.  140E+00  -9.  39AE-01  1  31 26*00  9  2346+01 
438E-01  1.  43AE+01  -8.  0846-01  1.  1336+00  1.  132E+02 
B13E-01  -2  3076*01  -1  202E+00  -3  A77E-01  7  2AOE+01 
83AE-01  -4. 1136-01  -1  777E-01  -3  8496-01  9  9426+01 
234E-01  —3.  3A9E-01  -4  A93E-01  -2.  349E-01  1  0046+02 
3A1E+00  -A.  020E-01  -A.  1076-01  2.  1326-01  1.0036*02 
1AAE+0O  -1.  AAAE+OO  -A.  177E-01  4  93AE-01  9  8736*01 
OOOE+OO  2  817E+0O  -3.  374E-0I  A  718E-01  1  0336+02 
88AE-01  -7  443E+00  -3  019E-01  A  402E-01  9.  2736+01 
1406-01  1  39A6+01  -3  331E-01  3  1A3E-01  1  1416+02 
3396*00  -2  343E+01  -7  474E-01  2.  A39E-01  7  341E+01 
2316-01  -2  494E-01  -3  9A8E-02  -2. 24AE-01  10096 +02 
2936-01  -1  397E-01  -1  3896-01  -A.  1436-02  1  0076+02 
2456-01  -3.  3336-01  -2  2106-01  8.  233E-02  1  004E+02 
847E-01  -1  3026*00  -2  4146-01  2  9496-01  9.9376*01 
3126-01  2  9936*00  -2  3006-01  3  8936-01  1.0336+02 
717E-01  -7  3016*00  -2  281E-01  3.1996-01  9.3026*01 
897E-01  1  3946*01  -2  A31E-01  2  2926-01  1  1436*02 
3896-01  -2  3386*01  -3  4AOE-01  -2.  2486-02  7  4406*01 
71AE-01  -1.7176-01  -1  A87F-02  -9  7896-02  99386+01 
4036-01  -1  4936-02  -9  8316-02  -1  0146-01  1.0006*02 
3206-02  -1.  A92E-0 1  -9  2896-02  A.  3946-02  9.  9886*01 
82A6-02  -1 .  3A9E+00  -1  0A4E-01  l  493E-01  9.8326+01 
3886-03  3.0476*00  -1  119E-01  1  871E-01  1  029E+02 
89AE-02  -7  3026+00  -1  191E-01  1  927E-01  9.2226*01 
8186-02  1  3886+01  -1  44A6-01  7  3116-02  1.  13AE+02 
931E-01  -2.  3316+01  -I  9A36-0I  4  2386-02  7  4436*01 
GA2E-01  -1  4826-01  -1  1A1E-02  -3.  0A3E-O2  1  0046*02 
0076-01  8  A24E-02  -2  9306-02  7.  0976-04  1  003E+02 
3096-01  -7  8986-02  -4  213E-02  7.3286-03  1.0046*02 
81A6-01  -1  298E+00  -3  IS3E-02  7.  032E-02  9.  913E+01 
925E-01  3  071E+00  -3  A27E-02  9  IA46-02  1  0346*02 
0386-01  -7  493E+00  -A  3836-02  A  3236-02  9  2806+01 
3236-01  1  384E+01  -7  4376-02  8.2996-02  1  141E  +  02 
1466-01  -2  34AE+0I  -8  84A6-02  3  3A76-02  7  4A7E+01 
7616-02  -1  3476-01  1  841E-03  3  70A6-03  9  993E+01 
8946-03  1  4226-01  A  39AE-04  -2  1406-02  1  0036+02 
83AE-02  -7  3A4E-02  -3  891E-03  A.  307E-03  l  0006*02 
4306-02  -1  207E+00  -A.  28AE-03  1  977E-02  9  8886+01 
1936-01  3  0246*00  -7  9A96-03  2  4706-02  I  0306+02 
AOAE-Ol  -7.  4AA6+00  -1  0716-02  l  2936-02  9  2346*01 
1836-01  1  3846*01  -1  4366-02  1  3036-02  1  I38E+02 
A77E-01  -2  339E+01  -l  7836-02  -2  1726-02  7  43AK+0I 


NODE  DISPLACEMENTS 
NO  U  V 

1  -1  712E-21  -3  402E-2I 

2  -3  420c -23  -2  938E-0P 

3  1  402F-21  -2  042E-03 

4  -3  33SE-29  -1  821E-02 

3  3  9619-22  -2  043E-O2 

6  -3  6B3E-2P  -2  140E-02 

7  -3  B09E-21  -3  377E-02 

0  1  076E-20  -3  937E-02 

9  l  .  242E-20  -2  234E-01 

10  2  B86E-20  1  033E-2O 

11  -8  018E-02  -1.744E-02 
12-7  329E-02  -l  913E-02 

13  -7  480c -02  -1  691E-02 

14  -7  489E-02  -2  S77E-02 
13  -7  094E-02  -1  324E-02 

16  -7  371E-02  -6  460E-02 

17  -6  310E-02  -1  609E-09 

18  -7  243E-02  -1  316E-01 

19  -6  832E-20  -3  621E-21 

20  -1  4039-0?  -2  407E-03 

21  -1  938E-02  -1  073E-02 

22  -1  08Oc-O2  -1  393E-0? 

23  -1  669E-02  -2  306E-02 

24  -1  3329-02  -1  014E-02 
23  -1  045c-02  -3  467E-02 

26  -6  344E-03  1  911E-02 

27  -1  960C-02  -1  027E-01 

28  6.  100E-20  -7  436E-21 

29  -7  137E-03  -2  426E-03 

30  -1  3929-02  -4  909E-03 

31  -1  426E-02  -8  873E-03 

32  -1  319E-02  -1  740E-0? 

33  -1  1329-02  -4  193E-03 

34  -9  6906-03  -4  107E-02 

33-1  374E-02  2  873E-0? 

36  -2  0229-02  -1  002E-01 

37  -3  341E-20  3904E-21 

38  -3  3019  03  -1  686E-03 

39  -6  497E-03  -2.  7/2E-03 

40  -8  0729-03  -4  8276-03 

41  -7  716E-03  -1  296E-0? 

42  -4  7909-03  2  OO7E-03 

43  -7  3329-03  -3.  333E-02 

44  -7  937E  03  3.  421E-02 

43  -1  0419-02  -9.  143E-02 

46  3  787t-20  -9  127E-21 

47  -6  9909-03  -8  074E-04 

48  -2  433E-03  -l  319E-03 

49  -2  889E-03  -2  600E-03 

30  -3  032E-03  -9.  483E-03 

31  -3  1429  03  3  483E-03 

32  -3.  2629  03  -3  197E-02 

S3  -3  129E  03  3.  771E-02 

34  -7  643E-03  -9  04 IE- OP 

33  -2  480c- 20  3  384E-21 

36  -1  7S6E-03  -9  093E-04 

37  -2  123E-03  -3.  434E-04 

38  -2.  770E-03  -9.  3«7E-04 

39  -2  98SE-03  -7.  761E-03 

60  -3  132E-03  7  741E-03 

61  -3  3629-03  -2  982E-0P 

62  -3  947t-03  3  927E-02 

63  -4  876E-03  -8  768E-0P 

64  1  238E-20  -6  993E-21 

63  2  761E-04  -3  724E-04 

66  -1  7729-04  -7.  613E-0S 

67  -4  064E-04  -3  947E-04 

68  -3  3639-04  -6  367E-03 

69  -6  388E  04  8  644E-03 

70  -9  473E-04  -2  874E-02 

71  -1  2376-03  4  032E-02 

72  -1  440E-03  -8.  709E-02 

73  -3  131C-20  2  783E-21 

74  3  676E-20  -7.  744E-04 

73  -3  0709-20  1  S12C-04 

76  3  384C-20  -2  667E-04 

77  -8  342C-21  -6  169E-03 

78  -4  420E-20  8  B63E-03 

79  11306-19  -2841E-02 

80  -2  3209-19  4  074E-02 

81  -2.  263E-19  -8.  734E-02 


